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Suppose that a quantum circuit with K elementary gates is known for a unitary matrix U, and 
assume that (7™ is a scalar matrix for some positive integer m. We show that a function of U can be 
realized on a quantum computer with at most 0{mK + logm) elementary gates. The functions 
of U are realized by a generic quantum circuit, which has a particularly simple structure. Among 
other results, we obtain efflcient circuits for the fractional Fourier transform. 
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Let U he & unitary matrix, U E Z^(2"). Suppose that a 
fast quantum algorithm is known for U, which is given 
by a factorization of the form 



U = UiU2---Uk, 



(1) 



where the unitary matrices Ui are realized by controUed- 
not gates or by single qubit gates We are interested 
in the following question: 

Are there efRcient quantum algorithms for 
unitary matrices, which are functions ofU? 

The question is puzzling, because the knowledge of the 
factorization (|l|) of U does not seem to be of much help 
in finding similar factorizations for, say, V — C/^/"^. The 
purpose of this letter is to give an answer to the above 
question for a wide range of unitary matrices U. 

Our solution to this problem is based on a generic cir- 
cuit which implements arbitrary functions of U, assum- 
ing that ?7™ is a scalar matrix for some positive integer 
m. If m is small, then our method provides an efficient 
quantum circuit for V. 

Notations. We denote by U (m) the group of unitary 
m X m matrices, by 1 the identity matrix, and by C the 
field of complex numbers. 



I. PRELIMINARIES 

We recall some standard material on matrix functions, 
see H, Q for more details. Let [/ be a unitary matrix. 
The spectral theorem states that U is unitarily equivalent 
to a diagonal matrix D, that is, U — TDT'^ for some 
unitary matrix T. The elements on the diagonal of 
D — diag(Ai, . . . , A2") are the eigenvalues of U . 
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Let / be any function of complex scalars such that its 
domain contains the eigenvalues Ai, 1 < i < 2". The 
matrix function f{U) is then defined by 

/(C/)=Tdiag(/(Al),...,/(A2.))^^ 

where T denotes the diagonalizing matrix of U , as above. 

Notice that any two scalar functions / and g, which 
take the same values on the spectrum of U , yield the 
same matrix value f{U) — g{U). In particular, one can 
find an interpolation polynomial g, which takes the same 
values as / on the eigenvalues Xi . It is possible to assume 
that the degree of g is smaller than the degree of the 
minimal polynomial of U. In other words, V = f{U) can 
be expressed by a linear combination of integral powers 
of the matrix U, 



(2) 



1=0 
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where m is the degree of the minimal polynomial of the 
matrix U, and G C for i — 0, . . . ,m ~ 1. In order for 
V to be unitary, it is necessary and sufficient that the 
function / maps the eigenvalues Xi of U to elements on 
the unit circle. 

Remark. There exist several different definitions for 
matrix functions. The relationshop between these defi- 
nitions is discussed in detail in 1^. We have chosen the 
most general definition that allows to express the func- 
tion values by polynomials. 



II. THE GENERIC CIRCUIT 

Let [/ be a unitary 2" x 2" matrix with minimal polyno- 
mial of degree m. We assume that an efficient quantum 
circuit is known for U. How can we go about implement- 
ing the linear combination (j^)? We will use an ancillary 
system of /i quantum bits, where fj, is chosen such that 
2^^~^ < TO < 2^ holds. This will allow us to create the 
linear combination by manipulating somewhat larger ma- 
trices, which on input |0) ® jV') G C^^ ® C^" produce the 
state |0)(g)l/|V'). 
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FIG. 1: A quantum circuit realizing the block diagonal matrix 
A = diag(l,^7,C/^...,^72"-l). 



FIG. 2: Generic circuit realizing a linear combination V. The 
case fi — 2 is shown. 



We first bring the ancillary system into a superposition 
of the first m computational base states, such that an 
input state |0) (8) \^) G ® is mapped to the state 



m — 1 



— y 



(3) 



i=0 



This can be done by acting with a 2^ x 2** unitary ma- 
trix B on the ancillary system, where the first column 
of B is of the form . . . , 1, 0, . . . , 0)*. EfRcient 

implementations of B exist. 

Notice that there exists an efficient implementation of 
the block diagonal matrix A ~ diag(l, U,U^, . . . , ?7^^~^). 
Indeed, A can be composed of the matrices C/^'', < 
77 < /z, conditioned on the fi ancillae bits. The result- 
ing implementation is shown in Fig. ^ The state is 
transformed by this circuit into the state 



E 

i=0 



(4) 



In the next step, we let a 2'^ x 2^ matrix M act on 
the ancillae bits. We choose M such that the state (jj) is 
mapped to 



(5) 



It turns out that M can be realized by a unitary matrix, 
assuming that the minimal polynomial of U is of the form 
X™ — T, T G C. This will be explained in some detail in 
the next section. 

We apply the inverse A'^ of the block diagonal ma- 
trix A. This transforms the state (||) to 



m 



|fc)«)F|V') 



(6) 



fe=0 



We can clean up the ancillae bits by applying the 2^ x 
matrix B^. This yields then the output state 



|O)0F|^) = \Q)®J{U)\^). 



(7) 



The steps from the input state |0)(8)|V-') to the final output 
state |0)(8)F are illustrated in Fig. for the case fi ~ 2. 



The following theorem gives an upper bound on the 
complexity of the method. We use the number of ele- 
mentary gates (that is, the number of single qubit gates 
and controUed-not gates) as a measure of complexity. 

Theorem 1 Let [/ &e a 2" x 2" unitary matrix with mini- 
mal polynomial x™' ^ t. t Cz C Suppose that there exists 
a quantum algorithm for U using K elementary gates. 
Then a unitary matrix V = f(U) can be realized with at 
most 0{mK + m? logm) elementary operations. 

Proof. A matrix acting on /i e 0(log m) qubits can be re- 
alized with at most O(m^logTO) elementary operations, 
cf. [Q. Therefore, the matrices B,B\ and M can be 
realized with a total of at most 0(3m^ logm) operations. 

If K operations are needed to implement U, then at 
most UK operations are needed to implement Ai(J7), 
the operation U controlled by a single qubit. The reason 
is that a doubly controlled NOT gate can be implemented 
with 14 elementary gates and a controlled single qubit 
gate can be implemented with six or fewer elementary 
gates 0. 

We observe that 2'' — 1 copies of Ai ([/) suffice to imple- 
ment A. Indeed, we certainly can implement Ai(t/^ ) by 
a sequence of 2'^ circuits Ai{U). This bold implementa- 
tion yields the estimate for A. Typically, we will be able 
to find much more efficient implementations. Anyway, we 
can conclude that A and A'^ can both be implemented by 
at most 14(2'' — 1)K S 0{14:mK) operations. Combining 
our counts yields the result. □ 



III. UNITARITY OF THE MATRIX M 

It remains to show that the state (^ can be transformed 
into the state (^ by acting with a unitary matrix M on 
the system of n ancillae qubits. This is the crucial step 
in the previously described method. 

Let U he a unitary matrix with a minimal polynomial 
of degree m. A unitary matrix V ~ f(U) can then be 
represented by a linear combination 



(8) 
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We will motivate the construction of the matrix M by ex- 
amining in some detail the resulting linear combinations 
of the matrices U^V . From dsh, we obtain 



i+k 



(9) 



Suppose that the minimal polynomial of U is of the form 
'm{x) = cc™ — q{x), with g{x) — J27Lo^ 9i^^ ■ "^^^ right 
hand side of (|^) can be reduced to a polynomial in U of 
degree less than m using the relation J7™ = g{U): 



m — 1 



The coefficients j3ki are explicitly given by 

where P denotes the companion matrix of m{x), that is, 



/ 1 
1 



P 





50 91 92 



1 

5m- 1 / 

The 2'^ X 2'" matrix M is defined by 

C 



M = 



1 



where C = {fiki)k,i=o,...,m-i, and 1 is a (2^— m) x (2^— rn) 
identity matrix. Under the assumptions of Theorem ^, 
it turns out that the matrix M is unitary. Before prov- 
ing this claim, let us formally check that the matrix M 
transforms the state (^) into the state (^). If we apply 
the matrix M to the ancillary system, then we obtain 
from (0) the state 



m— 1 



m— 1 



I I Ij _l_ -4 I I l! J_ 

^ m— 1 m — 1 

^ fc=o i=0 

-. m — 1 

= Y.\k)®i]^v\i,) 



which coincides with (H), as claimed. 

Lemma 2 Let U he a unitary matrix with minimal pol 
nomial m{x) = x™ — r. Let V be a matrix satisfying 
If V is unitary, then M is unitary. 



Proof. It suffices to show that the matrix C is unitary. 
Notice that the assumption on the minimal polynomial 
m{x) implies that C is of the form 



C 



ai 

"0 



an 



\ 



\ rai TOL2 ■ ■ ■ Tam~i ao J 



that is, C is obtained from a circulant matrix by mul- 
tiplying every entry below the diagonal by r. In other 
words, we have 



C=([r]. 



i>j imodm 



0,...,m — 1 



where [rjiyj — t if i > j, and [T]iyj ~ 1 otherwise. 

Note that the inner product of row a with row b of 
matrix C is the same as the inner product of row a + 1 
with row 6-1-1. Thus, to prove the unitarity of C, it 
suffices to show that 



Sa,0 



(row a row 



0) = r Q,-_nQ 

3=0 



(10) 



holds, where 6a,o denotes the Kronecker delta and the 
indices of a are understood modulo m. 
Consider the equation 



E 

V i=0 




(11) 



The right hand side can be simplified to a polynomial in 
U of degree less than m using the identity r [/™ = 1. The 
coefficient of C/" in (|l^) is exactly the right hand side of 
equation ([lO|). Since the minimal polynomial of U is of 
degree m, it follows that the matrices JJ", J7^, . . . , 
are linearly independent. Thus, comparing coefficients 
on both sides of equation ([l^) shows (^0|). Hence the 
rows of C are pairwise orthogonal and of unit norm. □ 

A Simple Example. Let Fn be the discrete Fourier 
transform matrix 



2^"/2(g^p(_2^j H/2"))fe,fco„ 



with = —1. Recall that the Cooley-Tukey decomposi- 
tion yields a fast quantum algorithm, which implements 
Fn with 0{n^) elementary operations. The minimal poly- 
nomial of Fn is x^ — 1 if n > 3. Thus, any unitary ma- 
trix V, which is a function of Fn, can be realized with 
0{n^) operations. 

For instance, if n > 3, then the fractional power F^, 
a; S R, can be expressed by 

P,^ = ao{x)I + ai{x)Fn + a2{x)Fl + a^{x)F^, 

where the coefficients ai(x) are given by (cf. |7j): 

ao(x) = ^(1 + e") cos a;, ai[x) ~ ^(1 — ie") sinx, 
a2{x) = 1 + e*^) cosx, Q;3(a;) ~ \{-~f — ie*^) sina;. 
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In this case, is reahzed by the circuit in Fig. g with 
U — Fn and M = {aj-i{x))i,j=o^,,_3. The circuit can be 
implemented with 0{n^) operations. 



IV. LIMITATIONS 

The previous sections showed that a unitary matrix /([/) 
can be realized by a linear combination of the powers C/*, 
< i < TO, if the minimal polynomial m(x) of U is of 
the form — t, re C. One might wonder whether the 
restriction to minimal polynomials of this form is really 
necessary. The next lemma explains why we had this 
limitation: 



Proof. By assumption, U"^ = rl for some t G C. This 
means that the minimal polynomial m(x) of U divides 
the polynomial — t, that is, x'" — r = TO(a;)TO2(a;) for 
some TO2(x) e C[x]. 

We may assume without loss of generality that the 
function / is defined at all roots of x™ — t. Indeed, we 
can replace / by an interpolation polynomial g satisfying 
/([/) = g{U) if this is necessary. 

Choose any unitary matrix A G J7(2") with minimal 
polynomial m2{x). The minimal polynomial of the block 
diagonal matrix Ua — diag(C/, A) is x'" — t, the least 
common multiple of the polynomials m{x) and m2{x). 
Express /{Ua) by powers of the block diagonal matrix 
Ua: 



Lemma 3 Let U be a unitary matrix with minimal poly- 
nomial m{x) — — g{x), degg{x) < to. IJ g(x) is not 
a constant, then the matrix M is in general not unitary. 

Proof. Suppose that g{x) — X^ilo^ Si^^'- We may choose 
for instance V — J7™ = 9{U). Then the norm of first 
row in M is greater than 1. Indeed, we can calculate 
this norm to be \go\^ + \gi\'^ + • • • + |gm-ip. However, 
|(?oP = 1, because go is a product of eigenvalues of U. 
By assumption, there is another nonzero coefficient gi, 
which proves the result. □ 



V. EXTENSIONS 

We describe in this section one possibility to extend our 
approach to a larger class of unitary matrices U. We 
assumed so far that f{U) is realized by a linear com- 
bination (||) of linearly independent matrices J7'. The 
exponents were restricted to the range < i < m, where 
TO is degree of the minimal polynomial of U. We can cir- 
cumvent the problem indicated in the previous section by 
allowing m to be larger than the degree of the minimal 
polynomial. 

Theorem 4 Let U E Z^(2") be a unitary matrix such 
that J7™ is a scalar matrix for some positive integer to. 
Suppose that there exists a quantum circuit which im- 
plements U with K elementary gates. Then a unitary 
matrix V = f{U) can be realized with 0{mK -\-m^ logm) 
elementary operations. 



f{UA) = diag(/(t/), fiA)) = J2 «.diag(C/^ A^). (12) 



The approach detailed in Section [II yields a unitary ma- 
trix M to realize this linear combination. On the other 
hand, we obtain from (Ha) the relation 



fiU) = a,U^ 



1=0 



by ignoring the auxiliary matrices A'^ , < i < m. It is 
clear that a circuit of the type shown in Fig. |^ with /i 
chosen such that 2^~^ < m < 2'' implements this linear 
combination of the matrices t/*, < i < to, provided we 
use the matrix M constructed above. □ 



VI. CONCLUSIONS 

Few methods are currently known that facilitate the en- 
gineering of quantum algorithms. Linear algebra allowed 
us to derive efficient quantum circuits for f{U), given 
an efficient quantum circuit for U, as long as C/™ is a 
scalar matrix for some small integer to. This method 
can be used in conjuction with the Fourier sampling tech- 
niques by Shor the eigenvalue estimation technique 
by Kitaev and the probability amplitude amplifica- 
tion method by Grover pO[ |, to design more elaborate 
quantum algorithms. 



[1] A. Barenco et ai, Phys. Rev. A 52, 3457 (1995). 
[2] W. Ferrar, Fimte Matrices (Oxford University Press, 
London, 1951). 

[3] F. Gantmacher, The Theory of Matrices, vol. I (Chelsea, 
1960). 

[4] R. Horn and C. Johnson, Topics in Matrix Analysis 

(Cambridge University Press, New York, 1991). 
[5] R. Rinehart, Amer. Math. Monthly 62, 395 (1955). 



[6] D. DiVincenzo, Proc. R. Soc. London A 454, 261 (1998). 

[7] B. Santhanam and J. McClellan, IEEE Trans. Signal Pro- 
cessing 44, 994 (1996). 

[8] P. W. Shor, in Proc. FOCS 94 (IEEE Computer Society 
Press, 1994), pp. 124-134. 

[9] A. Y. Kitaev, Russian Math. Surveys 52, 1191 (1997). 
[10] L. Grover, in Proc. STOC (ACM, New York, 1996). 



